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Abstract: In non-relativistic field theories, quantum fluctuations give rise to dissipa- 
tive behaviour even at zero temperature. Here we use holographic methods to explore 
the dissipative dynamics of massive particles coupled to quantum critical theories. We 
present analytic expressions for correlation functions and response functions. The be- 
haviour changes qualitatively as the dynamical exponent passes through z = 2. In 
particular, for z > 2, the long time dynamics of the particle is independent of its 
inertial mass. 
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1. Introduction 

At a quantum critical point, physics is invariant under the dynamical scahng symmetry 

x^Xx , t^XH (1.1) 

z is called the dynamical exponent. Such theories have been studied in great detail, 
starting with [1, 2]. For z = 1, the theory is typically Lorentz invariant, with the scaling 
symmetry promoted to the full conformal group. Theories with z ^ 1 arise naturally 
in a number of condensed matter contexts. 

A gravity dual for non-relativistic critical points was first proposed in [3]. The (d+l)- 
dimensional background metric is known as the Lifshitz geometry and takes the simple 
form 

ds'^^-r'^'de + ^ + r^d^ (1.2) 

where we have set the curvature of spacetime to unity. Several physical phenomena 
have been studied in this background, including mechanisms for linear DC conductivity 
and power-law AC conductivity [5] and fermions at finite density [6, 7, 8]^. 

^Intriguingly, the Lifshitz geometry appears to exhibit Ouroborosian stability properties, with 
divergent tidal forces in the far infra-red leading to a stringy instability [9] which is sometimes resolved 
[10] and sometimes not, the latter case resulting in a background AdS2 x geometry [11] which can 
itself be unstable to form a Lifshitz geometry [5] which has divergent tidal forces in the far infra-red 
leading to... 
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Despite enjoying both energy and momentum conservation, when z ^ 1 certain 
processes exhibit a form of dissipation, even at zero temperature. This is seen most 
clearly in the ballistic motion of objects — whether lumps of energy created from the 
fields or massive probe particles — both of which slow down, seemingly experiencing 
a friction force. There is no mystery here. Indeed, for relativistic z — 1 theories 
the constant speed of the centre of mass is not due to momentum conservation but 
instead can be traced to invariance under boosts. For -2 7^ 1, energy and momentum 
are conserved but drain into the soft infra-red modes of the theory. The higher the 
value of z, the greater the number of these modes^ and the greater the friction. 

Moreover, for such "dissipative" processes, there is sometimes a qualitative difference 
between theories with z <2 and theories with z > 2. Perhaps the simplest example is 
the following: take a massive particle in the metric (1.2) and throw it in a direction x 
parallel to the boundary. Of course, the particle will follow a timelike geodesic, falling 
into the infra-red, r — )■ 0, of the geometry. The question we would like to ask is: how 
far does the particle travel in the x direction? This is a simple exercise. One finds that 
for 2; < 2, the particle reaches a; — > 00 as t — )■ 00. But ioi z > 2, the particle only 
travels a finite distance. 

From the perspective of the boundary theory, this experiment corresponds to con- 
structing a lump of energy and giving it a kick. As the bulk particle falls towards 
r — )■ 0, the lump disperses, while the motion in the x direction tracks the centre of 
mass of the lump. For any z > 1, the lump slows down. But for z > 2 the centre of 
mass never gets further than some finite distance^. 

A zero temperature friction force also occurs for massive probe particles. In the 
context of string theory cmbcddings, these massive particles are modelled by strings 
attached to D-branes which sit at some fixed radial position in the background. The 
dynamics of such strings were first explored in AdS black hole backgrounds as a model of 
quarks ploughing through the quark-gluon plasma [4]. In Lifshitz spacetimes, analogous 
string dragging calculations reveal the existence of a friction force at zero temperature 
[5, 12, 13, 14]. Furthermore, in the regime of hnear response, there is again a difference 
between theories with z < 2 and z > 2. This manifests itself, for example, in the AC 

^This can be seen, for example, in the single particle density of states. A particle with dispersion 
relation ~ fc^ in d — 1 spatial dimensions has density of states p{E) ~ £;-i+('^-i)/^. 

^Although this statement holds true for any theory with a gravity dual, it does not appear to 
be a generic feature of quantum critical points. The counter-example is simply the non-relativistic 
Schrodinger equation, viewed as a classical field theory with z = 2. The centre of mass travels with 
constant velocity. 
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conductivity due to massive charge carriers which becomes non-analytic for 2; > 2 [5]. 
(We will review the calculation behind this in Section 3) . 

There is a simple dimensional argument for the crossover in behaviour that occurs 
for z = 2 [5]. The kinetic term for a particle with inertial mass m is 

S^^Jdt mX-X (1.3) 

The dimensions are inherited from the ambient critical fields, with scahng (1.1). This 
means [X] = — 1 and [t] — —z. The mass then has dimension [m] — 2 — z and the 
kinetic term on the worldhne of the particle is irrelevant for 2; > 2. The result is that 
the inertial mass, m, should play no role in the long-time dynamics of the particle when 
z>2. 

The purpose of this short note is to elaborate on this z — 2 crossover. Our main 
result, presented in the next section, is an analytic formula for the correlation function 
{X{t)X{0)) for a massive particle coupled to a quantum critical point. We will find 
that the long time behaviour is indeed independent of the inertial mass when z > 2. In 
Section 3, we revisit the calculation of the response function for the particle presented 
in [5]. We confirm the zero temperature fiuctuation-dissipation theorem computed 
through the bulk variables. We further describe the diffusion of the particle at finite 
temperature. 



2. Fluctuation 

Throughout this paper we model massive particles in the Lifshitz background as strings^, 
suspended from a D-brane that fills the spacetime hypersurface at radial coordinate 
r — Tb- The string hangs from this D-brane, and extends down into the IR r — )■ of the 
geometry (1.2). While the endpoint of the string on the D-brane represents the position 
of the point particle, the string itself is to be interpreted as the strongly-coupled critical 
fields sourced by the particle. 

In static gauge, the dynamics of the string are governed by the Nambu-Goto action 

S = ^— f ' drdt r'-\ 1 + r^x< ■ x' - -^'x ■ 'x (2.1) 

"^Lifshitz geometries in supergravity and, importantly, type II string theories have been constructed, 
starting with [5, 15, 16, 17, 18]. The Lifshitz geometry also appeared in a different context in [19] and 
discussed in [20]. 
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with Neumann boundary conditions imposed at the D-brane 

x'in) = (2.2) 

The classical ground state of the string is simply the trivial solution x — constant. The 
energy of this straight string is 

E ' ^' 



27ra' z 

E is the energy cost to create the string. For 2; 7^ 1, this is not the same thing as the 
inertial mass m in (1.3). Indeed, on dimensional grounds [E] — z and [m] — 2 — z so 
we expect 

m ~ (2.3) 

In Section 3, we will compute the response function and make this relation more precise. 

Our goal is to understand the quantum fluctuations of the end point of the string, 

X{t)=x{t,n) (2.4) 

We do so by computing the two-point Green's function for the x{t, r) fields, evaluated 
at Tfe, namely {X{t)X{Qi)). For small fluctuations, we hnearise the equations of motion 
about the average configuration x{t, r) — 0. The Nambu-Goto action reduces to 

const — ■ / drdt \ r^^^x' ■ x' tX ■ x 



Ana' J \ ^ 

with the resulting equation of motion 



dr \ dr J ^ dt 

Since the dynamics in each of the transverse dimensions decouple from one another, 
we will drop the vector notation and focus on just a single direction, x. Because the 
equation of motion is linear and invariant under time translations, we can Fourier 
decompose the most general solution as an integral over frequencies. 

x{t,r) = ^c.(r) [a{u)e-'^' + a(a;)^e*"*] (2.6) 

where the modes are given by the J and Y Bessel functions in the combination 

U^{r) ^ ^^=-^{ji,i {—) {—)) (2.7) 
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We have chosen a particular normahzation of these modes which will be explained 
shortly. The coefRcient C^, which fixes the relative weights of the J and Y Bessel 
functions, is determined by the requirement of Neumann boundary conditions (2.2) at 
the D-brane. It is 



Cuj — 



Ji_i 

z 2 



Yi_i 

z 1 



At this stage, we turn to the quantum theory. The momentum conjugate to a;(t, r) is 

7r(i,r) 



1 1 N 



and we impose the usual equal-time canonical commutation relations 

[x{t, r), 7r(i, r')] = i5(r - r'), [x{i, r),x{t, r')] = [7r(i, r), 7r(i, r')] = 0. 

The reason for the cleverly chosen normahzation of (2.6) and (2.7) now becomes appar- 
ent. Taking the inverse Fourier transform of (2.6), and using various integral properties 
of the Bessel functions, we find that the equal time commutation relations above imply 
that the creation and annihilation operators are correctly normalized^, 

[a{u), a\u')] = 5(u - u'), [a{uj), a(uj')] = [a^(uj), a^(u')] = 0. (2.8) 

Finally, we define the vacuum state to be 

a(a;)|0) = V a; 

This is the vacuum appropriate for observers that are stationary with respect to the 
coordinates in (1.2). 

We now have everything in place to compute the two-point function of the end point 
of the string (2.4). Substituting r — rf, into the mode expansion (2.6) and engaging in 
some appropriate Besselology, we find 



{X{t)X{0)) = 



27r 



e-"^' {X{uj)X{0)) 



with 



{X{uj)X{0)) ^8a'z 



( ^ 



J? . (- 

z 2 \zrl 



z 2 \zr 



(2.9) 



(2.10) 



^In fact, there is more to this calculation than we're letting on. To do things carefully, we first 
impose an IR cut-off which renders the spectrum discrete, as in [21]. This allows us to correctly 
normalize the modes, which then merge into a continuum as the cut-off is removed, resulting in (2.8). 
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It is instructive to look at the behaviour of the correlation function at low frequen- 
cies, for here we first see the crossover in behaviour at 2; = 2 that we described in the 
introduction. Prom the standard series expansions for the Bessel functions, we have 
that for u! <^ r^, 

(XHA-(0,> J ^^7!" (2.11) 

Constant factors in these expressions have been suppressed for clarity. As advertised, 
the energy E of the particle and, by (2.3), the inertial mass m, does not enter into the 
low frequency correlation for 2; > 2. 

In the special case of 2; = 1 (pure AdS), the Fourier integral (2.9) can be evaluated 
exactly to give 

For other values of 2;, no analytic expression for the inverse Fourier integral can be 
found. However, it is possible to obtain some estimate of the behaviour of the integrals 
at large times by assuming that the dominant contributions to the integrals come from 
the low frequency part of the integration range (2.11). Using these expressions, we find 
the long time behaviour 



4, .,2 2 



{X{t)X{0)) 



E^—z\t\^—z l<z<2, 
\t\^ z>2 



A marked transition occurs at 2; = 2. Here the two-point function grows linearly with 
t, which is its maximum rate of growth. The minimum growth occurs for 2; = 1 and 
2; — )■ 00, both of which exhibit logarithmic behaviour. Moreover, for 2; > 2, the long- 
time correlation of the particles has the peculiar behaviour of being independent of how 
heavy the particle is. 



3. Dissipation 

We now turn to discuss the dissipation of momentum as the particle moves through 
the quantum critical bath. We will continue to work at zero temperature. However, at 
the end of this section we will briefiy discuss the thermal contributions to the diffusion 
of the particle. 
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We start by computing the linear response of the particle to an external force F[u), 

{X{u)) = xM^M (3.1) 

The essence of this calculation can be found in [5] (see also [21]). A force F{uj) — Ee~^* 
is added to the endpoint of the string by the application of an oscillating electric field 
E on the D-brane. The string action receives an additional boundary term from the 
associated gauge field A^, 



S = dt At + A-x 



r=ri, 



Varying this action results in a new contribution to the boundary condition for the 
string endpoint. The Neumann boundary condition (2.2) is now replaced by 

-V.) = ^^- (3.2) 

'b 

As in the previous section, we look for solutions to the equation of motion (2.5). We 
again restrict to motion in just a single direction, parallel to F, now with x{t,r) = 
x(r)e~*'^*. Solutions are again constructed of Bessel functions, but this time our bound- 
ary conditions differ, both on the D-brane and, crucially, in the infra-red. Let us start 
with the infra-red. Here we impose ingoing boundary conditions, appropriate to the 
computation of the retarded Green's function x(a;). These select the first Hankel func- 
tion, H^^^ — J + lY , over the second, 

-('.'■) = ;i^^ri5(i^)^"'"' 

This solution is supported by a force F{(jj) which can be computed using (3.2). It is 

T?l \ 1 l-'2/2 t7(l) [ ^ \ 

F(uj) — cur, ' HW [ 

^ ^ 2TTa' -z-h \zrl) 

Using (3.1), we then have the response function 

It is again instructive to look at the low-frequency expansion of the response function. 
For u! <^r^, the expansion takes the form 

^^'^^ " m(ia;)2 + 7(-ia;)i+2A + ... ^^'^^ 
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with 



(2 - z) rl 



1 



.z-2 



1 
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1 r(i-i) 
(2.)2A r (1 + 1) 



The form of the response function (3.4) has a very natural interpretation, m is identified 
with the inertial mass (1.3). As we anticipated earher, we see that it is indeed related 
to the energy as (2.3). However, perhaps more surprising, it changes sign at z — 2. 

The second term in the denominator of (3.4) is the self-energy of the particle. For 
z > 2, the self-energy dominates over the inertial mass at low frequencies. This is the 
manifestation of the simple dimensional analysis argument given in the introduction. It 
is responsible for the non-analytic power-law behaviour in optical conductivity observed 
in [5]. Note that 7 also changes sign at 2; = 2, and it does so in such a way that the 
ratio 7/m is continuous at 2; = 2 and the response function suffers no pathologies. 

The Fluctuation-Dissipation Theorem 

The fluctuation-dissipation theorem relates the correlation function (2.10) to the imag- 
inary, dissipative part of the response function 



where nB{uj) — (e^^ — l)'^ is the Bose-Einstein distribution and captures the dissipation 
due to thermal noise. The "-|-1" in (3.5) is the contribution due to quantum noise and 
in the zero temperature limit, ^ — >■ 00, is all that survives. 

Of course, it is straightforward to prove the fluctuation-dissipation theorem (3.5) in 
quantum mechanics through spectral decomposition. Our purpose here is to verify this 
result from the bulk perspective. This has been previously studied in the AdS black 
hole background in [21, 22, 23] where the black hole horizon and the associated finite 
temperature play a key role in the connection. There has also been related work on 
understanding current noise in the holographic context [24]. Here we verify that the 
result also holds for quantum dissipation at zero temperature in Lifshitz backgrounds. 



{X{uj)X{0)) = 2 [riBiuj) + 1] Imxiuj) 



(3.5) 



Indeed, the result is straightforward. Using a standard Bessel function identity, the 
imaginary part of the full response function (3.3) is 



As expected, this is proportional to the correlation function (X(a;)X(0)) given in (2.10). 
We have 

confirming the zero temperature fluctuation-dissipation theorem for strings propagating 
in a Lifshitz background. 

Thermal Diffusion 

For completeness, we finish with a discussion of dissipation at finite temperature. In the 
present context, this means looking at the string correlation functions in the background 
of a Lifshitz black hole. There is now a rather extensive literature on constructing 
black holes in Lifshitz spacetimes. Many of these constructions are numerical, starting 
with [25], while analytic solutions often involve a running dilaton breaking the scale 
invariance. (See, for example, [26]). 

Here we give these difficulties short shrift. Following [5], we simply work in a fixed 
background, paying no heed to the equations of motion that it solves. The {d + 1)- 
dimensional metric we choose is 

ds^ = -r^'f(r)dt^ + —— + r^dx ■ dx 

The precise form of /(r) will not concern us. We require only that /(r) = has a 
single zero at the position of the horizon: /(r) ~ (r — rh) for r ~ r/j. The temperature 
of the boundary theory is the same as the Hawking temperature of the black hole, 

f,= 1. 

(z + d-l)rl 

A massive particle in the hot Lifshitz bath is once again modelled by a string sus- 
pended from a D-brane placed at r = r^. Our goal is to understand the thermal 
fluctuations of the endpoint of the string as it undergoes Brownian motion. This is 
a rather straightforward generahsation of the results in [21, 22, 23] for asymptotically 
AdS black holes. Over long time scales, the endpoint is expected to wander diffusively, 
with the characteristic expectation value, now taken in the canonical ensemble, given 

by 

{:{X{t)-X{0))^:) ^2{d-l)Dt (3.6) 

Here the :'s denote normal ordering. They remove the UV-divergent quantum fluctua- 
tions, leaving just the finite fluctuations due to the thermal bath. Our interest here is 
in computing the temperature dependence of the diffusion constant D. 
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Rather than compute the correlation function (3.6) directly, we will instead invoke 
the fluctuation-dissipation theorem, now in the form of the Einstein relation. The 
diffusion constant can be related to the response function defined in (3.1) by 

Our goal, therefore, is once again to compute x{^)i iiow at finite temperature. The steps 
are exactly the same as the calculation that we just performed at zero temperature, 
but with one exception: we can't find exact solutions to the equations. The linearized 
equation of motion for the string in the black hole background is now 

|:('-'^Vm|)+-^x = o. 

subject to the boundary condition (3.2) on the D-brane. In terms of the tortoise radial 
coordinate r* — J drf{r)~^r~^~^, the equation of motion takes a Schrodinger-like form 

where the effective potential is 

Although we are unable to solve these equations analytically, all that we require to 
determine the diffusion constant is the low frequency behaviour of the response function 
x(a;). Following [21], we calculate it by a patching procedure, splitting the spacetime 
into three regions: 

• Region A: r ~ r/j, -C a;^ 

• Region B: r ~ r/j, F » u;^ 

• Region C: r ^ Vh 

In region A, we can drop the potential V{r) from the equation of motion. We choose 
the causal, ingoing solution with xa — ae~^'^^* asymptotics. This is 

--W-(i-(JtI^"'<^-0^«(-^)) 

In region C, the emblackening factor /(r) is approximately unity, so the solution is the 
same superposition of Bessel functions (2.7) that we found in the Lifshitz geometry at 
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zero temperature. We are interested in the small uj behaviour of the solution, so we 
focus on the leading terms in its series expansion. 



xc = c,(l+0{ui'))+c,{^;f*\l + 0(u,-)) 
Finally, in region B, we drop the u term from the equation of motion. The solution is 

= h{l + Oi.-)) - b^rt'f -^.^^ (1 + Oiu-)) . 

At the UV end of region B, r S> r/i, we have /(r) ~ 1 and 

XB»6,(l + 0(u,^))-^jM^(l + 0(a.')), 

Prom this, we deduce that the solutions in regions B and C are consistent if 



6i = ci, &2 = -C2 ( — ) 



z + 2) 

2z) 



At the IR end of region B, r ~ r/j, we need only use the fact that /(r) has a single zero 
to find 

where A; is a constant of integration. Patching solutions in regions A and B requires 

bi + b2k = a, 62 = r- 

Combining all of these results, we find that at lowest order in a;, 

2 

xcin) = a(l + 0{u)), x'cin) = ^ (1 + OH). 

Thus we arrive at an expression for the response function at low frequencies 

, , 2TTia' 

'h 

Using the Einstein relation (3.7), we learn what we set out to find: the dependence of 
the diffusion constant on the temperature is given by 

D ~ /3^-\ (3.8) 
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Having struggled with us through our patching pain, the reader will be delighted to 
learn that this result can also be derived by some simple dimensional analysis. (A 
similar argument is given in [21]). If we assume that the mean free path of the particle 
depends only on the temperature, then L^fp ~ T'^^^. Meanwhile, if the relaxation 
time similarly depends only on temperature then r ~ 1/T. In time t, the particle 
undergoes N t/r collisions and, assuming a diffusive mode of transport, travels a 
distance AX — y/NL^aip- The upshot is 

T 

This reproduces (3.8). 

In line with the general theme of this note, we find that there is again a difference 
as we cross z — 2. For z < 2, the rate of diffusion decreases with temperature. For 
z > 2, we have the more usual situation where the rate increases with temperature. 
However, given the dimensional analysis argument above, this crossover does not seem 
to be related to the irrelevance of the inertial mass which occurs at zero temperature. 
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